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The measurement by the BES collaboration of J/tp "fpp decays indicates an enhancement at the 
pp threshold. In another experiment BES finds a peak in the invariant mass of tt mesons produced 
in the possibly related decay J/ip — > 'yn'^n'T]'. Using a semi-phenomenological potential model 
which describes all the NN scattering data, we show that the explanation of both effects may be 
given by a broad quasi-bound state in the spin and isospin singlet S wave. The structure of the 
observed peak is due to an interference of this quasi-bound state with a background amplitude and 
depends on the annihilation mechanism. 
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I. INTRODUCTION 

The search for exotic states in the NN systems has 
been pursued for a few decades, but significant results 
have only been obtained recently. An indication of such 
states below the NN threshold may be given by the scat- 
tering lengths for a given spin and isospin state. How- 
ever, in scattering experiments, it is difficult to assess a 
clear separation of quantum states. Measurements of the 
X-ray transitions in the antiproton hydrogen atom can 
select some partial waves if the fine structure of atomic 
levels is resolved. Such resolution has been achieved for 
the IS states [l] and partly for the 2P states [1]. One 
can also use formation experiment methods to reach spe- 
cific states. In this way, an enhancement close to the 
pp threshold, has been observed by the BES Collabora- 
tion Q in the radiative decay 

J/tp -fpp. (1) 

On the other hand, a clear threshold suppression is seen 
in the decay channel J/ip —f ir'^pp. To understand better 
the nature of these pp states, one has to look directly into 
the subthreshold energy region. This may be achieved in 
the antiproton-deuteron or the antiproton-helium reac- 
tions at zero or low energies. Such atomic experiments 
have been performed, although the fine structure resolu- 
tion has not been reached so far Another way to 
look below the threshold is the detection of NN decay 
products. Recently the reaction 

J/ip 77r+7r"?7' (2) 

has been studied by the BES Collaboration This 
reaction is attributed to an intermediate pp configu- 
ration in the J^'~'{pp) = state which corresponds to 



spin singlet 5— wave state. A peak in the invariant me- 
son mass is observed, interpreted as a new baryon state, 
named X(1835). The interpretation of the peak as a new 
X(1835) has been questioned by the Jiifich group Q- The 
latter view is supported by our calculations, but we sug- 
gest the origin of the BES finding to differ from the pos- 
sibilities presented in Ref. [7]. It is argued here that the 
peak is due to an interference of a quasi-bound, isospin 
0, NN state with a background amplitude. The same 
quasi-bound state was found in Ref. [8] to be responsible 
for the threshold enhancement in reaction ([T|). 

The purpose of the present work is to discuss the 
physics of TV TV states produced in these J/ip decays and 
relate it to atomic experiments. In reaction ([1]) only 
three pp final states are possible, as a consequence of 
the J^*^ conservation. These differ by the internal an- 
gular momenta and spins. Close to the pp threshold a 
distinctly different behavior of scattering amplitudes is 
expected in different states. A further selection of states 
is possible, but one has to rely on the analyzes of the 
elastic and inelastic NN scattering experiments. This 
has been studied in Ref. [1] within the Paris potential 
model [1, [l^, [ll|, , which is also used in the present 
work. 

The final pp states allowed by P and C conservation 
in the "fpp channel are specified in Table I. These are 
denoted as ^^~^^Lj or ^^+'^''^^+^Lj, where S,L and J 
are the spin, angular momentum and total momentum 
of the pair, respectively, while / denotes the isospin. A 
unified picture and a better specification in the radia- 
tive decays is achieved, semi-quantitatively, with an ef- 
fective three-gluon exchange model jj. This description 
indicates the final "/pp state to be dominated by the pp 
^^Sq partial wave. In this wave the Paris potential gen- 
erates a 52 MeV broad quasi-bound state at 4.8 MeV 
below threshold. This state is named iViVs(1870). A 
similar conclusion has been reached by the Jiilich group 
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TABLE I: The states of the low-energy pp pairs allowed in the 
J/tp — > 7p|) decays. The first column gives the decay modes 
to the specified internal states of the pp pair. The J^^ for 
the photon is f . The second column gives the J^^ for the 
internal pp system, the last column gives the relative angular 
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FIG. 1: The final state factor qf \Tif\ for the J/ip decays 
into 7 and pp. The latest version of the Paris model generates 
a quasi-bound state of F = 52 MeV and 4.8 MeV binding 
energy, and is the most consistent with the data. 



although the Bonn-Jiihch potential does not generate a 
bound state in the pp ^^Sq partial wave Q. 

Under the assumption that the tt''", tt" and ry' are pro- 
duced in relative S waves, the reaction ([2]), if attributed 
to an intermediate pp as suggested by the BES group, is 
even more restrictive than the reaction ^ . It allows only 
one intermediate state the pp ^Sq, which coincides with 
the previous findings. The presence of an intermediate 
pp state in reaction ^ is possible but not granted. We 
show below that a more consistent interpretation is ob- 
tained with the dominance of the NN{1870) state which 
is a mixture of pp and nh pairs. 

The content of this work is as follows. Sec. [TTl con- 
tains a description of the final state pp interaction and 
is included here for completeness. The subthreshold NN 
scattering amplitude, needed to describe reaction is 
defined in Sec. IIIII Sec. [IV] gives the equation to be 
solved to calculate the amplitude of the meson forma- 
tion through the intermediate NN interaction. The re- 
sults are presented and discussed in Sec. |V]together with 
some concluding remarks . 



For any multichannel system at low energies, described 
by an S wave /C-matrix, the transition amplitude from 
an initial channel z to a final channel / may be described 

by 



1 + iqfAff ' 



(3) 



where Aif is a transition length, Aff is the scattering 
length in the channel / and 9/ is the momentum in this 
channel [l^. In the following, channel / is understood 
to be the pp channel. Within the same formalism the 
scattering amplitude in channel / reads 



A 



ff 



1 -I- iqfAff 



(4) 



For S waves at low energies Aif,Aff are functions of 
and the main energy dependence of the amplitudes comes 
from the denominators in Eqs. ^ and With large 
values of Re Aff > one may expect a bound (quasi- 
bound) state. For large Re A// < a virtual-state is 
likely, but one cannot determine these properties with 
absolute certainty unless a method to extrapolate below 
the threshold exists. This is particularly true in the pp 
case where the absorptive part Im ^// is large because of 
the presence of many open annihilation channels. Since 
the final photon interactions are believed to be negligible, 
the energy dependence observed in the J/^* — > jpp decay 
rate reflects the energy dependence in qf \Tif\ . 

Practical calculations also indicate an energy depen- 
dence in Aff and the Watson approximation, i.e., the 
constant Aif is not applicable in a broader energy range. 
One needs to use Eq. ^ and a weakly energy dependent 
formation amplitude Aif 1/(1 + q'jrf) as explained in 
Ref.Q, where a best fit value = 0.55 fm was found. 
Figure [1] displays sizable model dependence of qf \Tif\^ 
for the ^Sq calculated for four versions of the Paris po- 
tential model d, [13, [nl, [l2|- These versions followed 
the increasing data basis which, for the most recent case, 
includes antineutron scattering and antiprotonic hydro- 
gen data. The threshold enhancement is attributed to a 
strong attraction in this partial wave. It does not prove 
the existence of a quasi-bound state, but such a state 
is indeed generated by the model in the ^^Sq wave [§]. 
There are additional arguments to support this result 
which follow from light p atoms. The absorptive am- 
plitudes can be extracted from the atomic level widths. 
With the data from Refs. 0, @ such an extraction was 
described in Refs. [l3| and [l^. The data allow to obtain 
only an isospin-spin average but as indicated in Fig. ^ 
the existence of a quasi-bound state is consistent with the 
atomic data. The increase of subthreshold absorption 
is also supported by the atomic level widths in heavy 
p atoms [15|. In addition to the broad S*— wave state 
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FIG. 3: (Color online) The real Re V{r, E) potential for A^iV 
interactions in the ^^So state. 



FIG. 2: The absorptive parts of spin-isospin averaged Np 
scattering amplitudes extracted from the atomic level widths 
in H, De, "^He and *He p [13]. Squares: S waves and cir- 
cles: P waves. The bottom scale indicates the energy below 
threshold. The curves, calculated with the Paris 2008 poten- 
tial, give the amplitudes separately: a„(p) denote the np or pp 
S-wave amplitudes, respectively and 6„(p) the corresponding 
P-wave amplitudes. The strong increase of absorption in the 
pp S wave is attributed mainly to the So state. 



the Paris potentials generate a narrow '^^Pi quasi-bound 
state which arises in Paris 08 and Paris 99 potentials. 
It gains some support from widths in the antiprotonic 
deuterium as indicated in Fig. O 

The procedure outlined above is based on a simple form 
of the low-energy final state wave function "^pp. At large 
distances and for small q/r, it becomes 



Tff e^p{iqfr) 



[1 - Aff/r] 



1 + iqfAff 



(5) 

The right side of this relation expands the wave up to q^ 
terms. At shorter distance the wave function is no more 
directly related to the scattering matrix and depends on 
details of the interaction in channel /. Integrated over 
an unknown transition potential Vi f it generates the for- 
mation amplitude Aif in the transition amplitude Tif. 
Eq. ^ was used with the Paris [8|] and Jiilich flfi| poten- 
tials. These potentials also generated the Aff. Formulas 
([3|) or ([5|) are useful above the threshold but cannot be 
simply extrapolated to the subthreshold region. The dif- 
ficulty is related to the momentum qf = ^j2^f^f;Ef^fi 
where /ijvw is the reduced mass. Above, the threshold 
Epffj is the kinetic energy in the CM system, below the 
threshold Ejqf^ is negative and qf becomes imaginary. 
The outgoing wave exTp(iqfr)/r becomes exp(— \qf\ r)/r. 



It damps strongly the interaction term in Eq. ^ and a 
more precise description is necessary. We now turn to 
this point. 
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FIG. 4: (Color online) The absorptive — Im V{r, E) potential 
for NN interactions in the ^^5*0 state. 



III. OFF-SHELL NN INTERACTIONS 

For further calculations one needs the off-shell exten- 
sion of the scattering amplitude in the energy as well as 
in the momentum variables. The most general extension 
for S waves is given by 
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In the next step, Eq. 



is used at negative energies 



f{k,E,k') 



2tt 



Mr, k)VNN{r, E)^+{r, E, k'^ dr, 

(6) 

where '^'^{r^E,k') is the full outgoing wave calculated 
with the regular free wave tpoi^jk') = sin{rk')/{rk'). 
In this equation the momentum k' is not related to 
the energy E. The Fourier-Bessel double transform of 
f{k,E,k') would generate a nonlocal f{r,E,r') matrix 
in the coordinate representation. So involved calcula- 
tions do not seem necessary as the experimental data are 
rather crude. We resort to a simpler procedure, standard 
in nuclear physics (for an application in the antiproton 
physics see Ref. [l^)- The subthreshold scattering am- 
plitudes are calculated in terms of T matrix defined in 
the coordinate representation by 



2tt 



Mr, k'{E)) 



with k'{E) = y^2^j^pfE. The T{r,E) is a local equiva- 
lent of the nonlocal T matrix in the sense that matrix el- 
ements in the S waves fulfill the relation /(fc, E, k'{E)) = 
J dr ipoir, k)T{r, E)ipo{r, k' {E)) valid in a narrow sub- 
threshold region where the last integral is convergent. 
The Vjsifi{r,E) is the recent Paris interaction potential 
which is used in the Schrodinger equation to calcu- 
late 'i>~^{r,E,k'{E)). The potentials used are plotted in 
Figs. [3] and |4] and the resulting scattering amplitude is 
given in Fig. [5] 
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FIG. 5: The real Re T(E) and imaginary Im T(E) parts for 
NN scattering amplitude in the ^^So state. 

To describe the intermediate NN state we need also 
the Fourier transform of T(r, E) 



E ~ Ej^ff. For positive energies this equation is not 
practical due to zeros in the denominator that occur in 
f{r,E) [see Eq. ©] at multiplicities of k' = ir/r. One 
could nevertheless use it for k' < n /rmax where r^ax is 
the distance at which the potential is cut-off. One has to 
set rmax ~ 2 fm if one wants to extend the calculations 
up to energies of « 20 MeV above the NN threshold. 
The normalization of f{r,E) by ipoir, k' (E)) in Eq. ^ 
insures the convergence of the integral [Eq. ([H])] below 
threshold. In order to see the predictions of our model 
above threshold, we shall replace in Eq. ([7]) ipo{r,k'{E)) 
by unity for > 0. 

The relevant on-shell 5'-wave scattering amplitudes 
given by 



T{E) = J drT{r,E), 



(9) 



T{k, E)^ dr T{r, E) 



sin(Kr) 



(8) 



are normalized to the corresponding scattering lengths 
at the threshold. The results for T{E) are plotted in 
Fig. O One can notice a resonant behavior of the imag- 
inary part. Roughly, the structure of this amplitude is 
dominated by the weakly bound state NNg{1870) in this 
wave. However, the location of the bound state given as a 
pole in the complex energy plane ( Re £^ = 1871.7 MeV, 
F = 52 MeV) corresponds neither to the maximum in the 
Im T{E) which occurs at £^ « 1840 MeV nor to the max- 
imum of TT~^TT~r]' invariant- mass distribution that occurs 
at 1835 MeV. The interpretation of the X(1835) turns 
out to be more involved. 

The potentials which generate this state are plotted 
in Figs. [3] and ID The real potential ( Fig. [3]) contains 
a very weak attractive tail, a repulsive barrier, a strong 
energy dependent attraction in the 0.5— 1.0 fm range and 
a repulsive core [ll|. These features, modified by the en- 
ergy dependent annihilation potential (Fig. 2]) and the 
proximity of the threshold, generate a rather untypical 
NN scattering matrix in the subthreshold energy region. 
The width of the ^^5*0 bound state indicates some en- 
ergy dependence. Moreover, the bound state form-factor 
displays strong enhancement in the subthreshold region 
which is a typical phenomenon of the subthreshold ex- 
trapolations. Altogether a strong enhancement of Im T 
is generated in the region well below the actual binding 
energy. As discussed above this effect finds support in 
the widths of the p-atom levels indicated in Fig. (5] 



IV. THE INTERMEDIATE pp STATES 

We assume that the photon in reaction ^ is emitted 
before the annihilation into mesons has taken place, as it 
happens in reaction ([T]). A specific model for that process 
was suggested in Ref. [1| but it will not be needed here. 
We assume however that the formation of the NN pair is 
described by a source function Fij and the annihilation 
by another function Ff^mes- In this way the effect of 
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the intermediate NN interactions can be described by 
an amplitude for the meson formation 



(10) 

where G{p,p' , Epf^) is the full Green's function for the 
intermediate NN system. The form assumed for the an- 
nihilation amplitude is 



Ff,„,es{p', Q) = (exp(-(Q - p')' t))) , (11) 

where the angular average over Q is indicated by the 
brackets. This choice is motivated by simple model con- 
siderations and the simplest possible assumption that the 
two TT mesons in reaction ([2]) are correlated to the /o (600) 
(also named a meson). The mass of the latter is as- 
sumed to be 500 MeV in our calculations. The relative 
momentum (2^ of the final r]' and a mesons is denoted 
by Q while the Gaussian profile comes from quark re- 
arrangement models of annihilation which operate Gaus- 
sian wave functions. 

The Green's function in Eq. pO|l may be expressed 
in terms of the free Green's function Go and the NN 
scattering amplitude T as 



G — Go + Go T Go- 



(12) 



Now, with the scattering amplitude defined by Eq. ([7]) 
and Eq. ([S]) one obtains 

T„„es - J dp dp' F,j{p)Go{p,Ej^^) [<5(p-p') 

+ T{\p-p'\,E^^) Go{p',E^j^)] Ff^mesip')- 

(13) 

The first term in Eq. (|13p corresponds to a background 
amplitude with a non-interacting NN pair. The sec- 
ond one describes intermediate state interactions. The 
Green's function is Go{p,E) = 47r/[(27r)^(gj - p^)] and 
the normalization is chosen such that T in Eq. has 
dimension of length. Let us notice that below the NN 
threshold both qj and Go are negative. Below the quasi- 
bound state T is attractive (negative) and the interfer- 
ence in Eq. (fT3l) becomes constructive. This effect ex- 
tends the peak structure to lower energies. We assume 
the formation amplitude to be described by 



1 + p-^r: 



2J1 



(14) 



with the range parameter — 0.55 fm determined before 
from the final interactions above the threshold [8]. The 
normalization is arbitrary. The angular integrations in 
Eq. (|13p generate an amplitude which depends only on 



IQI; that is due to the momentum dependence of the half- 
off shell T matrix and to the absence of any preferred 
direction in the initial NN state. 

A semi- free parameter is related to the radius pa- 
rameter in the quark models for the nuclcon and mesons. 
The range of allowed r/ values is limited. The upper limit 
r/ « 0.55 fm is obtained assuming the r.m.s. radii of the 
quark densities to be equal to the electromagnetic radii 
(0.8 fm for baryons and 0.6 fm for mesons). A lower limit 
r/ « 0.25 fm is obtained with the radii used in NN inter- 
action models based on quark approaches [17 [and quark 
rearrangement models of NN annihilation lg|. These 
rely on r.m.s radii in the range 0.5 — 0.6 fm for baryons 
and 0.4 — 0.6 fm for mesons. 

The last factor needed in this calculation involves the 
four body phase space for J/-0 ^ 77r+7r~r7'. We follow 
Ref. [13 to find 

dLi{M^j;P^,P^+,P^-,P^,) 



{Mj ^ Sm) 
(27r)24M2 



dL3{SM;P-,T+,P'7T-,Pri')dSM, (15) 



where Mj is the mass of J/^ and dL^ is the invariant 
phase space for the three meson system of invariant mass 
squared Sm- The dL^ may be found in Ref. [l^l and 
it generates only a weak energy dependence. The full 
phase space is used but one finds a simple approximation 
dL4 ~ e/ {rrirj' + 2m-,^ + e)^ with e = ^/Sm — rar,' — 2m7r 
to work well in the whole region of interest. 

All together the spectral function representing the 
X(1835) is given by 



dL^/dSM- 



(16) 



V. RESULTS AND CONCLUDING REMARKS 

The best description of the BES data is obtained with 
rf « 0.4 — 0.5 fm, and the shape of X(1835) calculated 
in this way is given in Fig. [51 The intermediate state 
is the isospin state. The data are reproduced fairly 
well despite the fact that the bound state itself occurs 
at 1871.7 MeV, i.e., 4.8 MeV below the threshold. This 
shape is determined by the interference effect of the two 
terms in Eq. describing the decay process. Within 
the Paris potential model and within a broad range of 
semi- free ri,rf parameters one finds no peak structure 
with the intermediate pp state. This result is consistent 
with the observation that isospin 1 for the final mesons 
is not allowed and the decay NN{T = 1) — > 7r+7r~?7' is 
not permitted by the isospin conservation [iS]. 

Other contributions, possible improvements. Above 
the NN threshold our estimation of Xs from our model 
Eq. Uni) using Eq. ([HI) with Eq. Q where ^po{r,k'{E)) 
is replaced by unity [see our discussion in Sec. M just 
below Eq. ([5])] generate a minimum (see Fig. [SJ that is 
deeper than the minimum indicated by the data. Below 
we indicate several possible explanations. 
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FIG. 6: The spectral function X5 representing the X(1835) 
shape. Here the range parameter of the annihilation am- 
phtude [Eq. Hill 

is rj = 0.45 fm. This S-wave contribu- 
tion has been normalized to reproduce the data close to the 
X(1835) peak. The experimental points are from Ref. [^. 
Above the A^A'' threshold the calculation is performed replac- 
ing ipo{r,k'{E)) by unity in Eq. (01. 
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and Xp (Eq. ^ with Ep = 1900 MeV and Fp = 200 MeV) 
and their sum representing the X(1835) shape. Here The S- 
and P-wave contribution have been normalized to reproduced 
the data ral. 



• The J^^ conservation allows the NN pair in ^^5o 
state to decay into the /o(600) ?/ pair in a relative S 
wave. This case has been discussed so far. In addition, 
with the baryons in ^^Pi states another decay mode to 
the /o(600) f]' pair is possible. It requires the two final 
mesons to be in the relative P-wave state. In the recent 
Paris 08 @ as well as in the former Paris 99 potential 
a close to threshold resonance is generated in the related 
^^Pi state. With the energy Ep = 1872 MeV and width 
Fp = 20 MeV it may contribute a spike to the spectral 
distribution in Fig. [S) Since two different partial waves 
are involved in the final states such a decay produces no 
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(17) 



to be added to the main expression for Xs given in 
Eq. The Xp possibility is a speculative one and 

the relative strength Cp would be very hard to predict. 
Also, with the recent update of the Paris potential the 
position of P wave resonance is not generated at "the 
proper" position. 

• In a more complete study, outside the scope of the 
present work, one could extend our 5'-wave equations 
[Eqs. Jll), (HI) and HH)] to the P-wave case. Here we 
illustrate in Fig. [7] the possible effect of an effective P 
wave represented by a resonant term given by Eq. (jl7|) . 
It can be seen that such an effective resonance with 
Ep = 1900 MeV and Tp = 200 MeV can fill up the above 
threshold dip of Xs- 

• The off-shell extension in terms of Eq. ^ cannot be 
fully trusted and the procedure of Eq. ^ should be used. 

• The final state factor given by Eq. pT|) is perhaps too 
simple to be used above the NN threshold. In some 
decay models, an energy dependent phase factor Pmes 
is expected [la]. This would have very limited effect 
below the threshold since the loop integrals - over Go - 
generate real functions. However, above the threshold 
the loop integral becomes complex and the interference 
pattern seen in figure (6) might be changed. 

These effects go beyond the technique used in this pa- 
per. 

• To confirm experimentally a direct link between the pp 
system and the X(1835), authors of Ref. Q have sug- 
gested a search at the future GSI Facility for Antiproton 
and Ion Research (FAIR) project using the proton an- 
tiproton detector array (PANDA) in reactions such as 
pp —> 7r"'"7r~A" and X Tr~^Tr^ri'. Another possible re- 
action would be pp jXilSSb). It could be performed 
with the PAX apparatus [2l| with ~ 50 MeV polarized 
antiprotons on polarized protons at CERN antiproton de- 
celerator (AD) Ring. The shape of the X{1835) could be 
tested by the photon energy distribution. Of special value 
would be the comparison of two measurements obtained 
with the parallel and anti-parallel initial spin configura- 
tions. That could give an information on the mechanism 
of the X{1835) formation. In particular, it would check 
the simple model presented in Ref. Q, where the ini- 
tial state (in the J/Psi case the intermediate) of the pp 
system is the spin triplet, which, after the emission of a 
magnetic photon, turns into the final spin singlet. 

In summary: it is shown that the X(1835) structure 
can be generated by a conventional NN potential model. 
Such a structure stems from a broad and weakly bound 
state, the iViVs(1870) that exists in the "^o wave. The 
existence of a quasi-bound S'-wave state receives an addi- 
tional confirmation from the level widths of antiprotonic 
atoms. 
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